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I. INTRODUCTION 

The Random Phase Approximation (RPA) is one of the most popular microscopic approaches to 
describe collective excitations in interacting many body problems pi 0, 0, 0, H|- For example it is 
widely used in condensed matter physics UIUQ, mainly to calculate the excitation spectrum via the 
linear response function. However, RPA is equally common in nuclear physics and other systems of 
finite size like cold atoms in traps @, El 0> or metallic clusters 0, etc. In these various fields RPA 
has slightly different meaning. We here will use it in the most general sense, like it is mostly employed 
e.g. in nuclear physics, including exchange, as it appears in linearising the time dependent Hartree-Fock 
equations around equilibrium 0, . Besides its relative simplicity, both conceptually and for numerical 
applications, RPA also has very desirable properties: when it is evaluated in the Hartree-Fock (HF) 
basis, then the f-sum rule and conservation laws are fulfilled, even and in fact most importantly in 
cases with spontaneously broken symmetries. For example in nuclear physics, where HF always breaks 
translational invariance, the corresponding RPA restores the symmetry in producing a Goldstone mode 
at zero energy 0,0. The same happens when other symmetries, like rotation (deformation) and particle 
number symmetries (superfluidity) are broken. Unfortunately when going beyond the HF-RPA scheme, 
the situation with respect to symmetry restoration and conservation laws is much less clear. Though 
there exists the well-known concept proposed by Kadanoff and Baym of the ^-derivable functional Q, 
which formally maintains all these properties, in practice any approach beyond HF-RPA faces serious 
difficulties. This stems essentially from the fact that in this scheme going beyond RPA is equivalent to 
introduce an integral kernel depending on more than one energy variable corresponding to the s,u, and 
t-channels of relativistic field theory. This then complicates the numerical task enormously. 

On the other hand, in certain circumstances one has to go beyond RPA because this approach also 
has a number of inherent quite severe drawbacks. For example it is well known that RPA violates the 
Pauli principle (the famous "quasi boson approximation"). This can induce strong errors and even give 
rise to qualitatively wrong results, as for example in the multilevel pairing Hamiltonian with Kramer's 
degeneracy, which is very much used recently for the description of superconducting nano grains |10| . 

A related problem which immediately appears when giving a close look to RPA is its lack of consistency. 
For example in the bubble summation one calculates the correlation energy of an electron gas Q but 

in the RPA-equation occupation numbers njj. = {0, 1} are used which correspond to an uncorrelated 
(HF) ground state. It seems obvious that, at least for the occupation numbers, a correlated ground state 
should be used. The corresponding approach already will contain RPA correlations self-consistently. 
This has become increasingly popular in recent years, and it has been called renormalised RPA (r-RPA) 

ana, Eaj . Lately r-RPA was widely and successfully used in describing the double beta decay process 
flEl fig and the electronic properties of the metal clusters . 
In fact r-RPA is an approximation to a more general extension of the RPA which is called Self-Consistent 
RPA (SCRPA) [HI, or also Cluster Hartree-Fock (CHF) [l8|. 

We recently had quite remarkable success with SCRPA in a number of non trivial problems |Tot [T^. l20| . 
In this paper we will specifically investigate the properties of SCRPA with respect to conservation laws 
and restoration of symmetries, the fulfilcmcnt of which is, as already mentioned, one of the outstanding 
aspects of standard RPA theory. It will be shown that for a wide class of symmetries SCRPA can be 
formulated in such a way that it also maintains these properties of standard RPA. This is to be considered 
as a strong advantage over other extensions of RPA because SCRPA, though much more demanding 
numerically than standard RPA, leads to an (nonlinear) equation of the Schrodinger type which seems 
to be accessible for a numerical solution. The restoration of symmetries within SCRPA is directly linked 
to the consideration of the so-called " scattering" terms and we will discuss their significance in detail. 

Here we will demonstrate those properties of SCRPA in an exactly solvable three level Lipkin model 
pH I22] ] . which can also be interpreted as an interacting spinproblem with two sites. Recently this model 
was investigated within the multistep variational approach [23 and the r-RPA formalism [24[ • This model 
has the property that, as swon in Ref. [25j, once the two upper levels become degenerate, a continuously 
broken symmetry develops beyond a certain critical coupling which can be considered as a deformed 
state with broken rotational symmetry. We will show that SCRPA, as standard RPA, will develop a zero 
mode (Goldstone mode) beyond a critical coupling, thus exactly restoring the original invariance of the 
Hamiltonian. 

We also will see that SCRPA yields much improved solutions which we will compare with the ones 
obtained from an exact diagonalization. For the two particle case SCRPA even reproduces the exact 
results. In detail the paper is organized as follows. In Section 2 we give a short description of the SCRPA 
and emphasise the main differences with respect to r-RPA. We outline here and in the Appendix B the 
method to determine the one and two-body densities. In Section 3 we give some numerical examples 
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concerning especially the Goldstone mode. In Appendix A we describe a very efficient procedure to solve 
SCRPA as a system of coupled equations, based on r-RPA as a first step. In the last Section we draw 
our conclusions. 



II. THEORETICAL BACKGROUND 



Let us consider a many body system of fermions characterized by n + 1 single particle levels, labelled 
a = 0, 1, 2, n. We define the " quadrupole-like" basis operators of this model 

N 

M =l 

where cL, denotes the creation operator of a fermion on some a-th level. We suppose the levels degenerate, 
the total number of projections /j, being N = 2fl. The operators (|2.1(l satisfy simple commutation rules, 
namely 

[K a 0,Ky$] = Sp-fKaS — SasKj/3 . (2.2) 

We will consider the following general Hamiltonian built on these operators 

n n 

H = 22 e aK aa + ^af3 1 sK a pK 1 s ■ (2.3) 

a=0 a/37(5=0 

We have chosen in our numerical application the three-level Lipkin model, with n = 2, corresponding to 
the SU(3) algebra. This model has been widely used in order to test different many-body approximations 
pll I22I I231 124. 12H . In analysing this model we have used a particular form of the Hamiltonian (|2.3() , 
namely 

2 2 

# = ^ e Q if M - - ^(^o^co + K 0a K Qa ) . (2.4) 

a— a— 1 

It is a particular case of the Hamiltonian used in Ref. j2|| and it is analyzed in detail in Ref. in 
connection with a continuously broken symmetry which appears, as already mentioned, when e\ —62- 

We will analyze the SCRPA in order to obtain the excitation energies and we will compare the results 
with the r-RPA scheme as well as with the exact results, obtained by the diagonalization pro cedure. 

We recall that both SCRPA and r-RPA schemes involve several common steps, namely [2|| 

A) One finds the minimum of generalised mean field equations coupled to the RPA fluctuations. 

B) One determines the RPA eigenstates. 

C) The RPA amplitudes are used to compute one and two-body correlated densities. At this step r-RPA 
supposes a factorisation for the two-body densities, while SCRPA uses a fully correlated approach. 

The amplitudes are used to iterate the steps A, B and C until the convergency is achieved. 
In the following we will shortly describe this procedure. 



A. Mean field 



As a first step in both SCRPA and r-RPA schemes we set up the generalised mean field equations. To 
this purpose let us introduce fermion creation operators in a general single particle basis as follows 

n 

where the supposed real transformation matrix is orthogonal. The unity matrix for Ck a corresponds to a 
"spherical" minimum, while a non-diagonal transformation corresponds to a "deformed" one. According 
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to Ref. |2jj for the three-level Lipkin model the above real matrix can be written as a product of two 
rotations, i.e. 

cost sincj) 0\ 1 1 
Cka = | sirup cos0 cosV' simp 
1 / \ — simjj costp; 

cost/) sirup cosip sum/) simp 

—sirup cos0 cosip cos0 sin?/; | . (2.6) 
— siiT0 cos?/' 

The basic operators ()2.1(l are replaced in this new basis by similar combinations of the new operators 
(1231) 

N 

By using the inverse transformation of 1)2.5(1 the Hamiltonian ((2.3(1 takes the following form in the new 
basis 

n n 

H = 2^ FijAij + ^ GijkiAijAki , (2.8) 

ij — ijkl—0 

where the coefficients are given respectively by 

n 

«=0 
n 

Gijkl = VafS^sCiaCjpCk-fClS . (2.9) 

The expectation value of the Hamiltonian ((2.8(1 with the correlated vacuum |0) has the form 

n n 

n = (o\h\o) = ]T F iS {ij) + ]T G miW) , (2.10) 

ij— ijkl—0 

where we introduced one and two-body densities as follows 

(ij) = (o\Aj\o) , 

(ijkl) = (0\A i:j A kl \0) . (2.11) 
In order to derive the Mean Field (MF) equations let us now consider the following restricted functional 

W = H- y £ j E a (aa) CkaCka ■ (2-12) 

a k 

We below will show that the terms with non-diagonal densities (aj3) vanish because they contain gener- 
ators of the RPA basis, which can be inverted in terms of SCRPA annihilation phonons acting on the 
vacuum state. The extremum condition 

|£--0 P-13) 

leads to the following system of equations 

e a (nn)C na + y^^y^^njkl;a(3~f5)C j f 3 Ck~ f Ci5 = E n (nn)C na , (2-14) 
M Pit 

where we introduced as short-hand notation 

(njkl; a/3jS) = ^[(njkl)V a pjs + (jnkl)V,3 aj s + (kjnl)Vy0 a s + (ljkn)Vs(3 ia } . (2-15) 



5 



By inserting a formal summation S nfl and using the orthonormality condition for the MF amplitudes 
one obtains an eigenvalue problem 

y^g nm C ma = E a (aa)C na . (2-16) 

m 

We solve this system to determine the eigenvalues E a and eigenvectors C ma for a = 0, 1, N. In case 
we evaluate expectation values H2.10(l using the standard HF vacuum one obtains what we shall call the 
"static minimum". If one uses the SCRPA correlated vacuum the result is called "dynamic minimum" 
(coupled with RPA fluctuations). In both cases the system l|2.16() contains nonlinear relations, because 
the matrix H nm 

H nm = (nn) e^Cn^Cm^ + (njkl; nf3jS)C m uC 3 0C k7 Cis (2-17) 

H jkl /j.PfS 

is written in terms of the eigenvectors and therefore it should be solved iteratively. For the "dynamic 
minimum" the iterative process is more complex because it involves the correlated vacuum, i.e. RPA 
amplitudes. As an initial solution we can consider the standard HF problem, leading to a decoupled form 
for the two-body matrix elements 

(ijkl) w (ii)(kk)Si 3 5 M + (m)(1 - {kk)/N)6 u 6 jk . (2.18) 
With (ii) ee pf ] = {1,0}, one obtains directly from Eq. I|2.14|l the standard system of HF equations 

^ 1 H n 2 l C ma = E a C ma , (2.19) 

m 

where H n °m are the elements of the HF Hamiltonian, which we will give for a particular interaction 
considered in our numerical application. However, in general, it is not necessary to consider the (kk) in 
(|2.18|l as uncorrelated. With correlations included Eq. I|2.18[l will lead to the so-called renormalised RPA 
(r-RPA) (see below). 



B. SCRPA equations 

The procedure to derive the system of RPA equations is given in many textbooks, see e.g. Q, 
Here we will shortly recall the main steps in order to introduce the necessary notations. To this purpose 
we use the equation of motion technique. The SCRPA creation operators are defined by the following 
superposition of the basic pair generators 

Qt = ]T (X^SQl, - Y r »JQ mt ) . (2.20) 

m>i 

These operators are written in the MF basis as follows 

SQL = Nj /2 A mi , SQ mi = 5Q\ m = N m ] /2 A im . (2.21) 
and define the excited states 

W) = Qt\0). (2.22) 

The RPA annihilation operators define a vacuum state, i.e. 

Q„|0} = 0. (2.23) 

It is important to stress that we used the letters m and i not for labelling holes and particles, but only 
to remember the condition m > i as written in H2.20|) . In this way we include in the basis all possible 
particle- hole (ph) and also the so called "scattering" elements with particle-particle (pp) and hole- hole 
(hh) configurations. This type of amplitudes is absent in standard RPA and we will have to discuss their 
significance in the course of the present paper. It has been recently shown that the inclusion of scattering 
terms allows to fulfil the energy weighted sum rule in the renormalised RPA (r-RPA) |2jj and they may 
therefore be of importance. 
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The normalisation factor in (|2.21() is given by the following mean value on the selfconsistent vacuum 

(0|[j4j m , A n j]|0) = S mn Sij((ii) - (mm)) = 5 mn 5ijN m i ■ (2.24) 

We then evaluate the mean value of the double commutator between the Hamiltonian and the phonon 
operator l|2.20|l 28] . One obtains in this way the SCRPA system of equations 

A B \ (X v \ _ (X v \ 
-B* —A* I \ Y v I ~ u}v \ Y v I ' 

where the matrices are given by the well-known relations 

Ami,nj 

(0| \SQ mi , IH^QljW |0) 



B mi , nj = -<0| [SQ mi , [H, SQ nj }} |0) = -Amijn ■ (2.26) 

Additionally we have the generalised mean field equations l|2.13[) . which with 12.2311 can also be written 
as 

(0| [H, Q v ] |0) = (0\[H, Qt] |0) = (0\[H, SQl] |0) = . (2.27) 
The amplitudes satisfy usual normalisation and completeness conditions 

E, x v — V v \ — A 

y ] (X m iX n j — Y rni X n j^ = S m i t nj . (2.28) 
v 

For the general Hamiltonian l|2.8|l one obtains the following form of the SCRPA matrices 
a . . = N 1/2 N 1/2 A ■ ■ 



■^{SijF mn - 5 mn Fji)((ii) + (jj) - (mm) - (nn)) 



1 

2 l 

^2(G ma bn(iabj) + G a bmn(abij))s + y^^G ai jb(amnb) + G a b 3 i(abnm)) s 



^(G majb (ianb) + G albn (ambj)) s , 
ab 

y^ y {Ggbcn(abcm) + G abmd (abnc)) s 

abc 

o <W ^ (Ggbci (abcj) + G ab j C (abic))s , 



abc 

byni^nj Qmi.jn ) (2.29) 

where we used the following short-hand notation 

(G ijk i(abcd)) s = G ljk i(abcd) + G kHj (cdab) . (2.30) 

The r-RPA differs from the SCRPA procedure in the way to determine the two-body densities. By 
using the factorised ansatz (|2.18J) one obtains the r-RPA matrices 



A 



(°) = hg. .jy(o) _ S H (0) )(N 1/2 N~ 1/2 + N~ 1/2 N 1/2 ) 



(a - —r \ n 1/2 n 1/2 



B (0) = l (s . H (o)_ S . H {Q) )(N 1/2 N~ 1/2 -N~ 1/2 N 1/2 ) 

G m inj - -j-jG nim ,j J N^N^ 2 , (2-31) 



7 



where we used the following notations 



Gminj ^mmj ^njmi ; 



(2.32) 



and 



11 mn — r mn i / y 



Gmaan ~t~ [ Gaamn Gmaan 



(aa) 



(2.33) 



With the help of the MF equation (|2~T§|) one obtains that -ffmn is diagonal, i.e. 

n mn — Vmn-^n ■ 



(2.34) 



We should stress that the occupation numbers appearing in (|2.33l) are calculated, as will be explained 
later, with the correlated vacuum and therefore (|2.34l) is not equivalent to standard HF, even within the 
r-RPA scheme. 

For r-RPA amplitudes with respect to the initial basis A m i 



0,mi — mi Q.mi ' fO,mi — mi 0,mi J 



(2.35) 



one obtains expressions with no phase space factors concerning the indices m,i, i.e. there appear no 
occupation numbers 



1 



h 0,mi 



lo^ — (E m — Ei, „ , 
-1 

U 0) + (E m - E^ 



E 

nj 

E 



Gr, 



G,, 



1 

N' 



-,G, 



Ig 



jy^jimn I * Q. n j 



G r , 



Gn 



1 

N' 



:G r . 



mm] I Y n j 



jyGnimj J X Q n j 



(2.36) 



where we defined the new amplitudes 



X 



— N rniXn 



0,mi — J "mi' l Q,mi ' 1 O.mi 



Y 



(2.37) 



We mention here that the r-RPA amplitudes defined by Eqs. H2.35J1 an d l|2.37|l are very important in 
solving numerically the SCRPA system of equations using the two-step method described in the Appendix 
A. 



C. SCRPA densities 



In order to estimate the expectation values of the one-body density operators within the SCRPA 

^ = ^4^, i = 0,1,2, (2.38) 

we need a special procedure. In the following we propose a method which is more convergent than the 
standard expansion in terms of phonon operators l|2.20|l . To this purpose let us introduce the shell model 
basis, but in the MF representation and which do not necessarily correspond to the HF basis, as we see 
in f2~Tfi|l 

Imna) = KU* A w A 20 \MF) , < m + n 2 < N , (2.39) 

where \MF) corresponds to a Slater determinant built with operators l|2.5|) . Within this complete basis 
we can expand a general combination of one-body densities 

N 

Po'Pi 1 ^ 2 = E b n 1 n 2 (kokik2)\n 1 n2)(n 2 ni\ 

nini- 
JV 

= £ ^W*>)^M^4' . ( 2 - 4 °) 

nini-0 
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where feo, ki, ki are some given exponents. The connection between the b and c coefficients is given by 
iterating the resolution of the unity. This method is used in Ref. [3(| and [3l| ■ For us it is important to 
find the coefficients c nin2 ( fecial ^2) and we will use a more direct technique. To this purpose we need the 
expectation values of the operators l|2.4U[l on the correlated vacuum |0). They are calculated by using 
the inversion of (|2.20() . One obtains 

Ann = N^^X^Q* + Y^Ql) , (2.41) 

where the normalisation factor N m i is given by l|2.24Jl . For instance in the case of products of two 
operators one obtains the following relations for to > i,n > j 

(OIA^-IO) = N^N^J2 Y ^ X nj 
(0\A un A ]n \0) = N^N^^X^ 

(0|A„„A„,|0) = N^N'J/Y^X^Kj ■ (2.42) 

V 

We mention that these relations can be directly used to express some of the SCRPA matrix elements in 
terms of RPA amplitudes. One finally obtains from (|2.4U|I a nonlinear system of equations, determining 
the normalisation factors A^iojA^o- It is given in the Appendix B in terms of expansion coefficients 
defined by Eq. H2.40JI . These coefficients are fast decreasing with increasing n\ + n-x- One can see from 
the Table 2 that already the linear approximation, i.e. n\ + ri2 < 1, ensures an accuracy of the order 
1CP 2 . In this case one also directly obtains the one-body densities 

N . 

(0|p |0) = JV-<0|fc|0>-<0|pa|0) , (2.43) 

where 

ymn = J2 Y ™oYn - (2.44) 



(0|p m |0) 



Vn 



2 3 

1 + (i/ii + 2/22) + "^2/112/22 



to = 1.2 



Based on this result one may hope that also for a more realistic n + 1-level model one can expand the 
product of densities in the ph "shell-model" basis, i.e. 

\ ni n 2 ...n n ) = AC 1 1 „ / 2 2 ..„„^ 1 ^...^ S |0), < n x + n 2 + ... + n n < 1 . (2.45) 



III. NUMERICAL APPLICATION 
A. Hartree-Fock mean field 

In order to apply the SCRPA and r-RPA procedures it is convenient to start iterations using standard 
HF and RPA solutions. To this purpose let us first evaluate the expectation values of the one and 
two-body operators on the standard HF vacuum 

(a/3) = SapSaoN , 
(apyS) = 6 aO 5 5o 5 M N[l + 5 0o (N-l)} . (3.1) 

We stress that the two-body density has the same form as the factorised ansatz (|2.18|l . The expectation 
value of the Hamiltonian l|2.4l) then becomes 

OnnO N + G 00 ooN(N - 1) 

a 

= Ne[eocos 2 4> + eisin 2 0cos 2 V' + e2sin 2 </>sin 2 'i/> 

- Xsin 2 0cos 2 0] , (3.2) 



9 



where we introduced the following dimensionless notations 



e k V(N-l) 

efe = — , X= 

e e 



(3.3) 



The Hamiltonian l|2.4|) has two kinds of HF minima, namely a spherical minimum and a deformed one 

1) = 0, V = , x < ei - e , 

(3.4) 



6j — CO 

2) cos 20 = , 0> = 0, % > ei — e . 



X 

Moreover, our calculations have shown that for any MF minimum one obtains ip = 0, defined by the 
parametrisation (|2.t)|) . independent of which kind of vacuum (correlated or not) we use to estimate the 
expectation values. 

B. Standard RPA 

For the above mentioned minima we obtain that the standard RPA matrix elements have very simple 
expressions. By using the following short-hand notation for the index pairs 



10 



1, 20 -> 2, 21 



(3.5) 



one obtains uL 



An 



e(ei — eo)cos 20 + e-^sin 2 20 , 



A22 
A12 = An 



e(e,2 — eocos 2 — eisin 2 0) + e— xsin 2 20 



= , 



sin 



B\\ = -ex(cos 4 <; 
#22 = -excos 2 , 
Bis = B 2 i = . 

Therefore the A and B RPA matrices are diagonal. The RPA frequencies are easy to estimate 



A 



kk 



B 



and the amplitudes become 



1 

71 



kk 1 



Akk 



± 1 



1.2 



1/2 



(3.6) 
(3.7) 

(3.8) 



We fix the origin of the particle spectrum with eo = 0. Then for a spherical vacuum = the RPA 
energies are given by 



1-1^ 



1/2 



f =1,2, 



with the corresponding RPA amplitudes 



1 
71 



± 1 



1/2 



(3.9) 



(3.10) 



As it was shown in Ref. J25|, if the upper single particle levels are degenerate, i.e. Ae = £2 — e i = 0, 
for the values of the strength \ > e i, m the "deformed region", i.e. with 7^ given by HF minimum, 
one obtains a Goldstone mode. In this case by considering e\ = 1 one obtains for the excitation energies 



wi = eV2(x 2 - 1) , 
U) 2 = . 



(3.11) 
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Indeed, in Fig. 1 by dashed lines are plotted standard RPA excitation energies versus the strength 
parameter \, by considering the number of particles N = 10. For x < 1 (spherical region) the eigenvalues 
u>i,u)2 are degenerate, while for \ > 1 (deformed region, 0^0) the second mode L02 has a vanishing 
value, i.e. it corresponds to the Goldstone mode. In this case the amplitude of the non-vanishing mode 
are given by 



ff) * ( 3 X-VX ±1 \ X i =Yi=0 . ( 3.i 2) 



The corresponding amplitudes for the Goldstone mode are given by 

x? = y x 2 = o, x| = y 2 2 ^ 00. (3.13) 

A similar plot for N = 20 is given in Fig. 2. 

C. SCRPA for N=2 

It is possible to find the exact eigenstates of the three-level Lipkin Hamiltonian (|2.4|) by diagonalizing 
it in the following normalized basis 



\nin 2 ) = f N N ZnT )l K ™ HF) ' ^ 

In this paragraph we will pay special attention to the spherical case </> = 0, with N — 2, i.e. the two- 
particle case. The above diagonalisation basis contains 6 elements with n\ + n 2 < 2. The calculation 
shows that the non-vanishing components of the ground state are \nin 2 ) — |00), |20), |02). The first 
two excited states contain only 1 10} and 1 01) components, respectively. This fact suggests that we can 
consider the two component phonon, with components (10) — > 1, (20) — > 2. The corresponding SCRPA 
vacuum state has the following form 

|0) = foo + ^K 2 W + ^Kl ) \MF) . (3.15) 



2 iu 2 

By using the action of the annihilation operator on the vacuum (I2.23|) one obtains for the coefficients 

-1/2 



«o 



ai = a Q (YX- 1 ) n , 

a 2 = a (YX- x ) 22 . (3.16) 

The one-body densities and their products can be expressed in terms of the vacuum amplitudes as follows 

<0|p,|0) = 2a 2 k , (0|^|0) = 4al , k = 0, 1,2 , 
(0|p Pi|0) - <0|poA»|0> = (0|pip 2 |0) = . (3.17) 

The SCRPA matrix elements of the Hamiltonian l|2.4|) can be directly evaluated by computing double 
commutators expectation values with respect to the vacuum. It turns out that they are diagonal and 
therefore the two eigenstates are decoupled. This allows us to express in a simple way the coefficients 
entering the vaccuum state 



^O^ 1 )^!^ k = l,2. (3.18) 
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The SCRPA matrix elements have the following form 



An 


= e 


ei - 


- e -+ 


2zi - 
"* 1- 


VZ2 

A \ 


A22 


= e 


62 - 


-e 4 


2z 2 - 
■X !_ 


Yz{ 

z 2 . 


Au 


= A 


21 


, 






Bn 




1 


+ *? 


+ Z\ z 2 






£X 


1 - 


z\ 


? 


B22 




1 


+ 4 


+ Z 2 Zx 








1 - 


z 2 
"2 




B12 


= B21 = 


. 







(3.19) 

The matrix elements of the Hamiltonian with respect to the above vacuum give the ground state energy 
as follows 



£o = 2e 



e 1 zf + e 2 z\ - x(zi 



Z2) 



1 



4 + 4 



(3.20) 



The SCRPA equations lead to the following nonlinear system of equations for zi, Z2 

X 



[z{ - 2z\ + {z\ ~ Zl ) Z2 + 1] + z\ - Zl 



X 

2e 2 



, 

[4 - 2z\ + {z\ - z 2 ) zi + l]+4-%2 = ■ 



(3.21) 



We solved this system by using the Newton procedure with a precision of 10~ 8 . The results for both 
eigenvalues and amplitudes coincide with the exact values using the diagonalisation procedure, indepen- 
dent of the used coupling constant x an d single particle energies. In Table 1 we give the coefficients 
dfe, k = 0, 1, 2, the ground state energy and the excitation energies cj&, k = 1, 2 for % = 5, eg — 0, e\ 
I . e 2 = 2, e = 1 MeV. The two lines correspond both to the exact and SCRPA values and they fully 
coincide. Therefore it is possible to construct the SCRPA ground state and excitations, which actually 
coincide with the exact solution. The inclusion of the third component (21) in the structure of the vacuum 
state for N = 2 gives a vanishing coefficient and therefore it is not necessary to consider it in this case. 



Table 1 



The vacuum coefficients ak,k = 0,1,2 in the ground state energy and the excitation energies 

ujk, k = 1, 2 at x = 5, e„ = n MeV, n = 0, 1, 2 for the exact and SCRPA values. 





ao 


01 


02 


Eo 




L02 


exact 
SCRPA 


0.772576 
0.772576 


0.496911 
0.496911 


0.395228 
0.395228 


-5.773794 
-5.773794 


6.773794 
6.773794 


7.773794 
7.773794 



For N > 2 the equation (|2.23(l leads to a system where the number of conditions exceeds the number 
of coefficients. In order to find the SCRPA ground state by satisfying the condition (|2.23|l it is necessary 
to enlarge the basis, by considering e.g. for N — 4 two particle-two hole(2p-2h) excitations together with 
the (21) component. Recently similar conclusions were obtained for the many level pairing model in Ref. 

As a matter of fact it is easy to convince one self that SCRPA is able to give the exact solution of any 
two-body problem. The ground state of a general two particle system can be written as 



10) 



1+2 E 

Zp 1 h 1 p 2 h 2 a p 1 a h 1 o,p 2 ah 2 \MF) , (3.22) 



where the 2p-2h coefficients obey the following symmetry relations 

z p 1 h 1 p 2 h2 Zp 1 h 2 p 2 h 1 — z p 2 h 1 p 1 h 2 z p 2 h 2 p 1 h 1 ■ (3.23) 
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Applying the vacuum condition (|2.23() on this state one obtains 



- Y vh + E z phP'h-X;r h , ala h \MF) = , (3.24) 



ph \ p'h' 



Y ph = ^2 z php'h'X^,* h , for all v j (3.25) 

p'h' 

which implies that 

ZpWh , = [YX~% hplh , . (3.26) 

Since we have now expressed the exact ground state in terms of the SCRPA amplitudes, we also can 
deduce that the SCRPA excitation energies are the exact ones because the RPA destruction operator can 
be written as Q v = |0)(^|, where \v) are the exact state vectors of the two particle excited states. On 
the other hand Q\ — |^)(0| fulfils exactly the SCRPA equation of motion 1)2.25(1 . We therefore see that 
SCRPA solves a general two particle problem exactly as long as one can define a Fermi surface. This, 
however, should in principle always be possible in putting the two particle system in an external box 
potential with or without periodic boundary conditions and solve the SCRPA equations in the limit of 
an infinitely large box. 



D. SCRPA in the spherical region 



Let us first discuss the SCRPA results in the spherical region, i.e. the region where the generalised 
mean field equation (|2.14|) has only the trivial solution = 0. In comparison with standard HF this 
region is strongly extended. The content of the spherical region depends on the particle number. We 
have seen that for N = 2 the spherical region covers all values of Xi since we obtained the exact solution. 
For N = 10 or N = 20 the spherical region is typically extended by a factor of two. This comes from the 
selfconsistent coupling of the quantal fluctuations to the mean field. 

Let us now consider the definite example eo = 0, e\ = 1, en — 2 for N — 20. In Fig. 3 we show by 
dashed lines the SCRPA results for the excitation energies, compared with the exact ones (solid lines) 
and to standard RPA (dot-dashes). We see that SCRPA strongly improves over standard RPA and in 
fact first and second excited states are excellently reproduced up to x-values of about \ ~ 1-2. The third 
state has no analogue in standard RPA and must therefore be attributed to the scattering configuration 
(21). The SCRPA solution for the the third eigenvalue definitely seems to approximate the fifth exact 
eigenvalue in the range < x < 1.0. The fifth eigenvalue has mainly a 2p-2h structure and it seems 
natural that the inclusion of scattering terms allows to reproduce such states, since e.g. a p a p > in a boson 
expansion picture is represented by a 2p-2h configuration, i.e. 

4<v ~Y. B U B p'h , (3-27) 

h 

where B p h are ideal boson operators (see |5|, Ch. 9). The fact that SCRPA, via the scattering terms, is 
able to describe states of the 2p-2h type is very satisfying and in a way astonishing. Indeed the single 
particle energy differences which correspond to the a p a p / configuration are e p — e p < and in our example 
this gives ea — ei = 1. One should have thought that for \ — > the third eigenvalue goes to one. In fact 
it can be shown that for x strictly zero, the third eigenvalue is one. However, this is a singular point and 
for any finite value of x the third SCRPA root jumps to the values shown in Fig. 3. A closer inspection 
reveals that this latter feature is entirely due to the genuine two-body correlations contained in SCRPA, 
because in the r-RPA the third eigenvalue is completely modified and indeed stays close to e% — e\ = 1. 
We repeated the calculation for Ae = e-i — e% — 0.001. In Fig. 4 for N — 10 and Fig. 5 for TV = 20 we see 
that the scenario concerning the third eigenvalue stays practically unchanged and, therefore, this seems 
to be a quite robust feature. In Fig. 6 we show the r-RPA solution for the Ae = 0.001 case. Indeed there, 
in the spherical region, a very low-lying third eigenvalue appears which is of the order Ae and thus not 
distinguishable from the abscissa on the scale of the figure. We therefore see that the genuine two-body 
correlations have dramatic effect in order to restore the situation for the third eigenvalue. 
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It is also interesting to investigate the influence of the scattering terms on the other SCRPA eigenvalues. 
We therefore performed the calculation switching off the (21)-component in the SCRPA equations for 
N = 10 case. In Table 2 we give the results. It is seen that the influence of the scattering terms on the 
first and second SCRPA eigenvalues is very weak. 

Table 2 

The eigenvalues for the two-dimensional ui?\ k — 1,2 and three dimensional uj^\ k — 1, 2, 3 versions 

of the SCRPA versus the strength \ (first column) in the spherical region. The particle number is 
N — 10 and e n = n MeV, n — 0, 1, 2. In the columns 2-4 are given the exact solutions u>k , k = 1,2, 5. 
In the last columns are given ground state energies for two-dimensional and three-dimensional SCRPA 

and exact values, respectively. 



X 




L0 2 


UJ 5 


u> x 




ui 1 




w 3 


^0 


^0 




0.900 





745 


1 


758 


2.888 





774 


1 


908 





778 


1 


910 


3 


102 


-0.407 


-0.409 


-0 


354 


0.920 





735 


1 


752 


2.887 





767 


1 


905 





772 


1 


907 


3 


106 


-0.427 


-0.429 


-0 


371 


0.940 





726 


1 


745 


2.886 





761 


1 


902 





766 


1 


904 


3 


111 


-0.447 


-0.450 


-0 


388 


0.960 





716 


1 


739 


2.886 





755 


1 


899 





760 


1 


901 


3 


116 


-0.468 


-0.470 


-0 


406 


0.980 





706 


1 


733 


2.885 





749 


1 


896 





754 


1 


898 


3 


122 


-0.490 


-0.492 


-0 


423 


1.000 





696 


1 


727 


2.885 





744 


1 


893 





749 


1 


895 


3 


128 


-0.512 


-0.514 


-0 


442 


1.020 





686 


1 


722 


2.885 





738 


1 


889 





744 


1 


891 


3 


135 


-0.534 


-0.536 


-0 


461 


1.040 





677 


1 


717 


2.886 





733 


1 


886 





739 


1 


888 


3 


143 


-0.557 


-0.559 


-0 


480 


1.060 





667 


1 


712 


2.886 





728 


1 


883 





734 


1 


885 


3 


151 


-0.580 


-0.582 


-0 


500 


1.080 





657 


1 


707 


2.887 





723 


1 


880 





729 


1 


882 


3 


160 


-0.604 


-0.606 


-0 


520 


1.100 





647 


1 


703 


2.889 





719 


1 


877 





725 


1 


879 


3 


169 


-0.628 


-0.631 


-0 


540 



In Figs. 7 and 8 we show the SCRPA groundstate energies for N = 10 and N = 20, respectively as 
a function of the angle <j> representing the mean field transformation. We take Ae = 0.001. We see that 
for a large plateau of x-vahies the minimum stays at \ — 0- This region is more or less enhanced by a 
factor of two, compared with standard HF theory, where the = minimum ceases to exist at \ = 1. In 
SCRPA, after the = solution is finished, the system jumps in a discontinuous manner to the "new" 
deformed solution. This discontinuity is clearly seen in Figs. 9 (N — 10) and 10 (N = 20) where we 
show the ground state energy E ^ CRPA ^ (dashed line) and compare it with the exact values (solid line) 
and standard RPA (dot-dashed line). The full 3x3 SCRPA with the scattering terms was used. We see a 
strong improvement of SCRPA versus standard RPA in the transition region where standard HF switches 
from the spherical to deformed regime. Beyond this transition point SCRPA stays spherical and slightly 
overbinds whereas standard RPA goes deformed and slightly underbinds. From the point where also 
SCRPA jumps to the deformed solution, standard and SCRPA give almost identical correlation energies. 



E. SCRPA in the deformed region 



In the deformed region a particular situation arises in our model for Ae = 0, since, as already mentioned, 
a spontaneously broken symmetry occurs in this case. Here the standard HF-RPA exhibits its real strength 
because, as shown in H3.11|l . a zero mode appears which signifies that the broken symmetry is restored, 
i.e. the conservation laws are fulfilled Q, ||. We will now demonstrate that this property is conserved in 
SCRPA (with scattering terms). In this context we mention that the symmetry operator 

La = i{K 21 - K 12 ) = i [A 2Q - A 02 )sin</> + (A 21 - A 12 )cos0] , (3.28) 

which commutes with the Hamiltonian, can be identified with the z-component of the rotation operator. 
This operator contains scattering terms which are, however, also present in our RPA-ansatz and therefore 
we can expect that (|3.28ll is a particular RPA operator with a zero eigenvalue solution in SCRPA, very 
much in the same way as this is the case in standard RPA. The numerical verification of this desirable 
quality of SCRPA must however be undertaken with care. Indeed, a zero mode contains diverging 
amplitudes which, injected into the SCRPA matrix, may not lead to selfconsistency. The way to overcome 
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this difficulty is to start the calculation with a finite small value of Ae, i.e. with a slight explicit symmetry 
breaking, and then to diminish its value step by step. We, in this way, could verify with very high accuracy 
that the zero eigenvalue occurs in the deformed region for all values of the interaction strength x- This 
is shown in Fig. 4 by a solid line for N = 10 and Fig. 5 for N = 20. Here we considered the value 
Ae = 0.001, but we were able to reach the value Ae = 10~ 6 . 

In spite of this clear appearance of the Goldstone mode, which only shows up when scattering terms 
are included and the generalised mean field equation (|2.14|) is solved self consistently with the SCRPA 
equations, the spectrum in the deformed region needs some discussion and clarification. First of all one 
notices a discontinuity when passing from the spherical to the deformed region. This feature of deformed 
SCRPA has been already noticed earlier with other models: the two level Lipkin model [28| and the two 
level pairing model [29). It is reminiscent of a first order phase transition which, however, is absent in 
the exact solution. 

Furthermore the second SCRPA eigenvalue seems too low compared with the exact solution and also, 
astonishingly, with respect to the standard RPA solution, as seen in Figs. 1,2. In fact standard RPA 
works surprisingly well in the deformed region, a fact which seems to be a constancy in all the models 
we have investigated so far [2^.l29|. 

Third, in the deformed region it does not seem easy to identify the third SCRPA eigenvalue with any of 
the exact solutions. Therefore, in spite of satisfying the Goldstone theorem, the interpretation of SCRPA 
seems less clear in the deformed region than it appears in the spherical one. However, the situation may 
be less unsatisfactory than it appears at first sight. 

As usual with a continuously broken symmetry also in the present model a clear rotational band 
structure is revealed, as can be seen by inspecting Figs. 1,2. The exact solution, found by a diagonalisation 
procedure, has a definite angular momentum projection Lq. Moreover the expectation value of the L\ 
operator has integer values, namely 



sJ{lJ) = J = 0,1,2,... (3.29) 

The first "rotational band" J = 0, 1, 2, ... is built on top of the RPA excitation with a vanishing energy 
(Goldstone mode). The second "band-head" corresponds to the second RPA mode (dashed curve) and 
the corresponding members of this "excited band" have a different slope with respect to the "ground 
band" . On the same figures one can also see a third " rotational band" . 

As customary in RPA theory one also can evaluate the mass parameter of the rotational band within 
SCRPA. By a straightforward generalisation we obtain for the moment of inertia (see e.g. Ref. 0,0) 

M = 2L* Q {A-By 1 L , (3.30) 

where A, B are the SCRPA matrices and Lq denotes the part with m > i components of the momentum 
operator Ij3.28|) . which should be written in terms of normalised generators SQ^ defined by Eq. (j2.21|) . 
i.e. 



L 



o = i(o, N^ 2 sin<f>, N^cosA . (3.31) 



For the standard RPA case, by using the matrix elements given by (|3.6|l . and A^o = N, one obtains an 
analytical solution, namely 

A(x-l) 

M = —P -{ . 3.32 

exCc+1) 

In Fig. 11 we show by a solid line the above moment of inertia calculated with RPA for N = 20 and in 
Fig. 12 with SCRPA. The comparison with the "exact" mass parameter fitted from the exact spectrum 
is also given by dashed curves. From the dot-dashed lines, giving their ratio, we see that the SCRPA 
mass is slightly closer to the "exact" value than the corresponding RPA value. This is also manifest when 
comparing the rotational energies, as in Fig. 13. Here we plotted by dashes the SCRPA spectrum and 
by solid lines the "exact" values. 

One therefore sees that one can recover rather clearly the precise position of the lowest energy value 
and therefore the discontinuity, at least of the first excited state, in crossing the phase boundary, is much 
less pronounced, than the one which is seen in Figs. 4, 5. One may think that the same upward shift 
should also be observed for the second SCRPA eigenvalue (the band head of the second rotational band) 
bringing it in a close agreement with the exact solution. It then appears that one can use SCRPA with 
scattering terms in the same way as standard RPA with all usual properties conserved. Let us mention 
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that even in r-RPA the Goldstone mode appears when the scattering terms are considered. This is shown 
in Fig. 6. Note, however, that in the spherical region the third r-RPA eigenvalue also comes almost zero 
value. The reason for this was discussed above. 

F. The occupation numbers 

The reproduction of energies may not be the most sensitive indicator for the performance of a theory. 
Usually wavefunctions are much more demanding in this respect. We therefore calculated for Ae = 0.001 
the occupation numbers which reflect more directly the quality of the wave function, since they are not 
obtained from a minimisation procedure, as the energies. 

In Figs. 14, 15 we show the SCRPA results for (0|p |0), (0|pi|0) and (0|p 2 |0) according to Eq. |2~15|) . 
We see that in the spherical region the SCRPA densities (solid lines) are rather close to the exact values 
(dashes) for %-values up to slightly above one. This is, of course, much better than standard RPA (not 
shown) which breaks down, as the energies, well before x = 1, but the region 1 < \ < 2 remains not well 
reproduced. On the contrary in the deformed region the SCRPA results are excellent. It is somewhat 
unexpected to see that the occupation numbers of the two upper but degenerated levels are not equal 
but their averaged sum very well reproduce the exact occupancy of the two upper levels (dotted lines in 
Figs. 14, 15). Therefore the SCRPA wavefunction performs well in the deformed region. 

IV. CONCLUSIONS 

In this paper we studied various aspects of the SCRPA in the three level Lipkin model. With respect 
to the two level Lipkin model, the three level version has the advantage of allowing for a continuously 
broken symmetry on the mean field level with the appearance of a Goldstone mode. One major objective 
of the present work therefore was to investigate whether in SCRPA the very desirable properties of 
standard RPA concerning exact restoration of symmetry and fullfillment of conservation laws can also be 
maintained. We found out that this is indeed the case under the condition that the symmetry operator 
is fully included in the RPA operator as a particular solution. In our model this implies that the RPA 
operator contains in addition to the usual ph and hp components a^ah and a^a p also the so-called 

anomalous or scattering terms aj,a p < and a h ah>- 

With a correlated ground state the inclusion of such terms is perfectly possible since the occupation 
numbers are now different from one and/or zero and mathematically the situation becomes similar to 
the finite temperature RPA case, where also pp' and hh' components are allowed. However, it was found 
in this work that the inclusion of such scattering terms is not without problem. For example pp 1 {hh 1 ) 
configurations can imply very small energy differences e p — e p < w and therefore for small interaction very 
low-lying roots may be found in the spectrum which are absent in the exact solution. In our model it 
turned out that this is the case only in the approximate, so-called renormalised version of SCRPA. On the 
contrary, in the full SCRPA the eigenvalue which is linked to the inclusion of scattering terms corresponds 
to a 2p-2h state of the exact spectrum. This is true in the "spherical" region whereas in the symmetry 
broken phase the third eigenvalue does not seem to have any precise correspondent in the exact solution. 
In the symmetry unbroken phase the elimination of the anomalous terms has an insignificant influence 
on the reminder of the spectrum. In the symmetry broken phase this inclusion is crucial to produce the 
Goldstone mode at zero energy, otherwise this state is pushed up to quite high energy with no equivalent 
in the exact spectrum. We also calculated the moment of inertia of the rotational band which works in 
SCRPA (with scattering terms) in a way very analogous to standard RPA. Very good agreement with the 
exact solution is found. Therefore the present formulation of SCRPA allows to maintain all the formal 
and desirable properties of standard RPA. There is only one unpleasant feature which is that SCRPA 
solutions do not join smoothly when passing from the symmetry unbroken to the symmetry broken phase. 
Whether this can be improved upon remains to be seen in future work. A further very attractive feature 
of SCRPA is the fact that it solves any two body problem exactly. This is contrary to the usual where 
many body approximations deteriorate when decreasing the number of particles. 

Quite remarkably the SCRPA solution exists in the symmetry unbroken single particle basis well beyond 
the phase transition point given by the standard RPA and the excited state is well approximated in all 
this region. For N = 10 and N = 20 at a certain strength Xcrt ~ %Xcrt, where Xcrt is the value where 
standard HF becomes unstable, SCRPA also jumps to the deformed solution. A similar behaviour of 
SCRPA has already been found for other models. 
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Though the appearance of the Goldstone mode is a very desirable feature, since this signals that 
conservation laws and sum rules are fulfilled, it turned out that in the deformed region, contrary to 
the "spherical" one, the SCRPA only very slightly improves over standard RPA and both, SCRPA and 
standard RPA are quite close to the exact solution. Again this is a feature which also has been found 
in other models. It seems that the common wisdom that "deformation" sums already a lot of additional 
correlations in standard RPA is born out here. 

We therefore can summarize the studies of SCRPA in various models in the following way: in the 
symmetry unbroken regime SCRPA works excellently, yielding simultaneously very good ground state 
and excitation energies. One can work in the spherical basis beyond the usual transition point given by 
the instability of the standard HF solution. This "beyond" depends on the number of particles and the 
model. SCRPA solves a general two particle system exactly. 

In the deformed region SCRPA does not seem to improve very much over standard "deformed" RPA. 
However, all formal properties of standard RPA are maintained with SCRPA when the scattering terms 
are included. For example for nuclei SCRPA will give a Goldstone mode restoring translational invariance. 
The transition from "spherical" to "deformed" is discontinuous, which is an unpleasant feature. It is not 
evident that the inclusion of scattering components in the spherical region improves the results. It also 
should be mentioned that the appearance of the Goldstone mode in SCRPA is only guaranteed when 
the symmetry operator does not contain any diagonal components as ajet;, since those terms cannot 
be included in the RPA operator. Such a situation occurs for example in the case of particle number 
violation 
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Appendix A 
SCRPA in two steps 



The SCRPA system of equations (|2.25|) and the normalisation condition for amplitudes can be written 
in a matrix form respectively as follows 

sx = xn , 

X^X = 1 . (A.l) 
Here we introduced some short-hand notations 

s - (;«•)• *-(;$)■ 



Here the matrices a, b are defined by Eq. (|2.29|) . We used the amplitudes X relative to the initial 
generators A m i defined by Eq. I|2.35[) and also the SCRPA matrices defined by Eq. (|2.29() . We also 
introduced similar to 1)2. 37|) amplitudes X. N is the diagonal metric matrix NkSki, where we used for the 
basis indices the notation (mi) —* k 

(10) -> 1, (20) -> 2, (21) -> 3 . (A.3) 

SCRPA equations are written as 

f/SX = Xfl , (A.4) 
in terms of the SCRPA matrix relative to X amplitudes 

S=Af- 1 SAf- 1 . (A.5) 
From the normalisation condition one derives the inverse amplitude matrix 

X- 1 =lX f , (A.6) 

giving the following resolution of the unity 

Xf??=\= {\ fj . (A.7) 

Let us split the SCRPA matrix into two parts: the r-RPA matrix plus a fluctuation 

S = S + SS , (A.8) 

where SS = S — Sq. The r-RPA system of equations and the normalisation condition are given by 

SoXo = Xoflo , 
XqXq = 1 , 

X = JVXq . (A.9) 

We should stress on the important fact that our calculation is sclfconsistent because we use the SCRPA 
set of amplitudes in deriving the metric matrix Af in the r-RPA system of equations. We insert the unity 
operator written in terms of r-RPA amplitudes 

X ixl=l, (A. 10) 

in front of X amplitudes of the left hand side Eq. (jA.ip . By using r-RPA equations (| A. 9|> one gets 

(^o^o + SSX ) TXpC = Xfl . (A.ll) 
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We multiply this relation to the left with X Q = XX\. One finally obtains a standard system of RPA 
equations 

(n i + ss') x' = ix'n , (A. 12) 

with the normalisation condition for new amplitudes 

X' ] 1X' = 1 , (A. 13) 

where 

yi vt~v 

SS' = IXISSXoX . (A. 14) 

In particular one has 

5a' = xlSax + xlSby + y^Sbxo + ylSay , 

5b' = -x\5ay Q - x\5bx - yl5by - yl5ax . (A. 15) 

This is an RPA system of equations with r-RPA energies on the diagonal. It corresponds to a new 
representation of the SCRPA phonon 

Qt = £(«L-«o«) , (A. 16) 

a 

in terms of r-RPA phonons 

Qla = { x 0kAmi ~ Vak^-irn) ■ (A-17) 

The formulation of SCRPA, where the division by the norm matrix Af has disappeared, has the impor- 
tant formal aspect that one sees that the division by eventual very small numbers or even zero, contained 
in A/", in reality does not exists. This is also born out by our three component SCRPA solution treated 
in the main text. 
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Appendix B 
Operator expansion method for densities 



In order to find the coefficients of the expansion (|2.4l)[l it is possible to use the resolution of unity 
method [30l l3l| . We proceed in a more direct way by computing the matrix elements with respect the 
complete basis (|2.39|) . One obtains a lower triangular system of equations 



^2 X! ■^mxm 2 C nin 2 {k klk 2 ) = £ mi m 2 (k kik 2 ) , 
ni<mi n 2 <m 2 



where 

\Antn2 — Kf—1 /()\ A m 2 A m i A n i A n 2 A n 2 A n i A m i A m 2\(]\ 
Jvl m 1 m 2 — ■ /V mim 2 \ U I /1 02 ^01 A 10 "^20 A 02 A)l A W A 20 l U / ) 

The solution of this system is given by the following recurrent relation 

1 



c mi m 2 {kokik 2 ) 



C mim2 {kokik 2 ) - ^ M^l^Cn^ikok^) 



nin 2 <mim 2 



J v ^m\m 2 

Let us consider the HF ansatz of the vacuum, defined by 

The coefficients of the system (|B.l|l are given in a straightforward way 

C mi rn 2 (kokik 2 ) = (N -m 1 -m 2 ) ko m'l 1 m% 2 

un ina _ K r {N -mi -m 2 + m +ti 2 )! 

M mi m 2 ~ A/ mim2 M(mi _ ni ) !(TO2 _ n2)! 
N\mi\m 2 \ 



Nrn\m 2 — 



(N — mi — m 2 )\ 
coo(k kik 2 ) = N ko S kl0 6k 2 o ■ 
We give the solutions of this system for n\ + n 2 < 2 
T k , (N-l) k -N k 



P K « N* + 



1 



N 



-(A w Aoi + A 20 A 02 ) 



(N -2V + N"- 1 (N -2) 

2N(N-l) [ 
x {A 2 1Q A 2 01 + A 2 20 A 2 Q2 + 2Ai A 20 A 02 A i) 



2(N-l) k+1 



N 



Prn 
P0Pm 



1 1 

W A m0 A 0m + N{N _ l} 



N - 1 



N 



-A m oAo 



1 



1 



N(N- 1) 



N 2 (N - 1) 

^4-10^4-20^-02^-01 ) 



(2 fe 1 - 1 + jj)A 2 m0 A 2 0m + -Ai A 20 A 02 A i 



(^mO^Om + Ai A 2 qA 02 A i 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



where m = 1, 2. These expansions are very convergent. Indeed, this can be seen in Table 3, where there 
are given the expansion coefficients c ni „ 2 (k kik 2 ) with m + n 2 <2 for N = 10. 

In order to find the normalisation factors using the expansion H2.40JI up to n\ + n 2 < 2 we need 
the expectation values for the following products of operators 

(0\A mO A Om \0) = N m0 Y,YXo 

V 



yv yv yfi yfl , yv ytf yv 1 
1 ??7,0 1 mO 1 nO 1 nO "r 1 mO 1 mO 1 n0 1 nOJ 



(B.7) 
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For 7ii + ri2 < 2 one obtains the following system in terms of the expansion coefficients in Eq. I|2.40|l 

N w = cqo(IOO) - coo(010) + [cio(100) - c 10 (010)]iVioJ/u + [c m (100) - c i(010)]iV 20 y 22 

+ [c 20 (100) - c 20 (010)]jV 1 Vii + [c 02 (100) - c 02 (010)]iV 2 2 z 22 

+ [cu(100)-cu(010)]JVi JVao*ia , 
^20 - cqo(IOO) - c 00 (001) + [c 10 (100) - c 10 (001)]iVio2/ii + [cqi(IOO) - c 01 (001)]iV 20 y 22 

+ [C2 (100) - O2o(001)]JV? *n + [c 02 (100) - c 02 (001)]iV 2 V 22 

+ [c 11 (100)-c 11 (001)]7V 10 A r 2 o2i 2 . (B.8) 

Table 3 

The expansion coefficients c ni „ 2 (kokikz) in Eq. \2.4-Wjl with n\ + n 2 < 2 for N = 10. 



ni,n 2 


0,0 


0,1 


1,0 


0,2 


1,1 


2,0 
















1,0,0 


10.0000 


-0.1000 


-0.1000 


-0.0011 


-0.0022 


-0.0011 


0,1,0 


0.0000 


0.0000 


0.1000 


0.0000 


0.0011 


0.0011 


0,0,1 


0.0000 


0.1000 


0.0000 


0.0011 


0.0011 


0.0000 


2,0,0 


100.0000 


-1.9000 


-1.9000 


-0.0100 


-0.0200 


-0.0100 


0,2,0 


0.0000 


0.0000 


0.1000 


0.0000 


0.0011 


0.0122 


0,0,2 


0.0000 


0.1000 


0.0000 


0.0122 


0.0011 


0.0000 


1,1,0 


0.0000 


0.0000 


0.9000 


0.0000 


-0.0011 


-0.0011 


0,1,1 


0.0000 


0.0000 


0.0000 


0.0000 


0.0111 


0.0000 


1,0,1 


0.0000 


0.9000 


0.0000 


-0.0011 


-0.0011 


0.0000 
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FIG. 2: The same as in Fig. 1, but for N = 20. 




FIG. 3: SCRPA excitation energies versus the strength parameter x, for N = 20 and eo = 0, ei = 1, ei = 2 
(dashed lines). By solid lines are given the lowest exact eigenvalues and by dot-dashes the standard RPA energies. 
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FIG. 4: SCRPA excitation energies versus the strength parameter x, for N = 10, Ae = 0.001 MeV (full line). 
By dashes are given the lowest exact eigenvalues and by dot-dashes the standard RPA energies. After the phase 
transition point \ — 1 (standard RPA) and \ ~ 2.1 (SCRPA) a Goldstone mode at zero energy appears. 




FIG. 5: The same as in Fig. 4, but for N = 20. 
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FIG. 6: r-RPA excitation energies versus the strength parameter x, for N = 20, Ae = 0.001 MeV (solid lines). 
By dashes are given the lowest exact energies. The Goldstone mode appears at zero energy beyond x ~ 1-4. For 
X < 1.4 appears an unphysical mode of the order Ae which cannot be distinguished from the abscissa on the scale 
of the figure. 
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FIG. 7: The SCRPA expectation value of the Hamiltonian versus the angle cj>, for N = 10 and different values 
of the strength parameter x (from the top of the figure, \ — 0, 0.5, 5). By dashes are given the values for the 
spherical region and by solid lines for the deformed region. 
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FIG. 9: The expectation value of the Hamiltonian versus \ for N = 10. By the lower solid line are given 
exact values, by dashes the SCRPA energies, by dot-dashes standard RPA results and by the higher solid line 
HF values. Once the SCRPA jumps to the deformed solution at x ~ 2.1 SCRPA and standard RPA are almost 
indistinguishible on the scale of the graph. 




FIG. 10: The same as in Fig. 9, but for N = 20. 
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FIG. 11: The standard RPA inertial parameter versus x (solid line), the inertial parameter from the exact energy 
spectrum (dashed line) and their ratio (dot-dashed line) for N — 20. 
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FIG. 12: The SCRPA inertial parameter versus x (solid line), the inertial parameter from the exact energy 
spectrum (dashed line) and their ratio (dot-dashed line) for N — 20. 
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FIG. 13: The SCRPA rotational spectrum (dashed lines) and the exact energies (solid lines) for N = 20. The 
three levels correspond to J — 1, 2, 3. 
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FIG. 14: The one body densities versus x for N = 10. By solid lines are given SCRPA values and by dashes 
the exact results. The dotted line represents the average of two occupation numbers of states 1 and 2 for 
X > Xcrt w 2.1. 
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FIG. 15: The same as in Fig. 14, but for N = 20. 



